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HAMILTONIAN RELATIVE EQUILIBRIA WITH CONTINUOUS 

ISOTROPY 


JAMES MONTALDI AND MIGUEL RODRIGUEZ-OLMOS 


Abstract. In symmetric Hamiltonian systems, relative equilibria usually arise in 
continuous families. The geometry of these families in the setting of free actions of 
the symmetry group is well-understood. Here we consider the question for non-free 
actions. Some results are already known in this direction, and we use the so-called 
bundle equations to provide a systematic treatment of this question which both 
consolidates the known results, extending the scope of the results to deal with non¬ 
compact symmetry groups, as well as producing new results. Specifically we address 
questions about the stability, persistence and bifurcations of these relative equilibria. 
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1. Introduction 

This article deals with a particular kind of integral curve of Hamiltonian flows on sym- 
plectic manifolds equivariant under the canonical action of a Lie group. They are usually 
known in the literature as relative equilibria, or steady state solutions, a terminology 
inherited from analytical mechanics. In this article we follow an approach to the qualita¬ 
tive study of some local properties of these solutions, based on the use of geometrically 
adapted tubular neighbourhoods. This approach is an evolution of methods previously 
proposed in the literature by several authors (see [23, 24, 26]). 

Let (VyUj) be a symplectic manifold, G a Lie group acting on V by symplectomor- 
phisms and h a smooth G-invariant function on V. Hamilton’s equations produce a 
vector field Xh on V, and due to the invariance of h and the symplectic structure under 
the G-action, the flow of this vector held sends group orbits to group orbits. A relative 
equilibrium is an integral curve of Xh that belongs to a single group orbit. 

The importance of relative equilibria is more apparent using orbit reduction. The 
flow of Xh naturally descends to a continuous flow on the quotient space VjG, and 
relative equilibria project to hxed points of this reduced how, therefore being (together 
with periodic orbits, which will not be addressed in this article) the primary object of 
qualitative local studies. If the action of G on P is free and proper the quotient space 
VjG is a manifold and the reduced how can be shown to be smooth and Hamiltonian 
with respect to some reduced Hamiltonian function and Poisson structure on the quotient 
space. This means that relative equilibria can be studied via their projected hxed points 
using many of the standard available techniques from differential geometry and critical 
point theory on this quotient space. 

For these reasons, this article deals with equivariant Hamiltonian hows for which the 
symmetry Lie group action has singularities. In this situation this geometric approach 
is no longer possible since, to start with, the quotient space VjG is no longer smooth. 
We focus on the problems of nonlinear stability of relative equilibria, as well as on the 
organization of relative equilibria in parametrized branches and their bifurcations. The 
study of these topics in the presence of singularities has been going on for a number of 
decades now, and a variety of methods have been used to attack them, from the use of 
singular reduction and conhnement arguments [15], to the use of suitable linearizations 
of the Hamiltonian held [6, 5, 2] passing by topological methods [10, 13, 21], to cite a 
few. Analogous questions for free actions are addressed in a number of papers, and in 
particular [19, 10, 20] 

The motivation for this article is twofold. On the one hand, we wish to obtain a 
convenient framework adapted to the study of the local dynamics of equivariant Hamil¬ 
tonian hows in the presence of singularities specihcally adapted to the unique geometric 
features of the problem, and powerful enough to serve as a unifying approach to deal 
with all aspects related to this topic. For that, we have adopted the “bundle equations” 
proposed in [23] and [24] and have built on this formalism in order to obtain a self 
contained geometric machinery especially adapted to deal with Lie group actions with 
singularities. On the other hand, we show how these ideas are implemented by proving 
within this framework some standard results in the theory, with improved hypotheses, 
as well as obtaining several new relevant ones. 
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Organization. This article is organized as follows: Section 2 collects some standard 
generalities on symmetric Hamiltonian systems and relative equilibria. Then the bundle 
equations, the main theoretical framework of this article, are introduced. We present a 
modification of the bundle equations of [23] and [24] which explicitly incorporates the 
existence of continuous stabilizers for the group action, resulting in the bundle equations 
with isotropy (2.10), (2.11) and (2.12). These are then used, in Proposition 2.6, to 
obtain a system of differential equations locally characterizing relative equilibria for a 
given equivariant Hamitonian flow. One of the main objectives of this article is to find 
solutions of these equations, and for that we prove two technical results (Lemmas 2.8 
and 2.9) that will apply in the subsequent persistence and bifurcation results. 

In Section 3 we use the bundle equations with isotropy in order to revisit the proof of 
well known stability results from [5, 15] and [12]. Our approach allows us to extend these 
to the statement in Theorem 3.6. Compared to [5, 15] we are able to avoid making use 
of the “orthogonality” choice of velocity, and compared to [12] we allow the possibility 
of non-compact momentum stabilizer G^. In Remark 3.7 we explicitly show with an 
example how our approach gives sharper stability results than those in [5] and [15]. 

In Section 4.1 we deal with the problem of persistence. That is, under which conditions 
a relative equilibrium persists to a continuous branch of relative equilibria. This differs 
with the problem of bifurcations, which considers as starting point a parametrized branch 
of relative equilibria and finds conditions for which at a particular value of the parameter, 
new relative equilibria exist arbitrarily close to the original branch. These new relative 
equilibria are typically organized in bifurcating branches. 

We obtain four main results in this direction: Theorems 4.3, 4.4, 4.5 and 4.6, based on 
different non-degeneracy or degeneracy hypotheses. Some of these results, or particular 
cases of them, had been obtained in the past using quite different approaches. Here 
we show how they all follow in a straightforward way from the bundle equations with 
isotropy and are therefore unified under a common framework. Many of these results 
assume the group of symmetries to be compact, whereas we allow non-compact symmetry 
groups in all our results, thereby also providing minor extensions to the known results. 

Finally in Section 5 we address the problem of finding bifurcations from parametri¬ 
zed branches of relative equilibria. Using the notion of branches of relative equilibria of 
the same symplectic type introduced in Definition 3.4, we are able to give, in Theorem 
5.1, sufficient conditions under which one such branch exhibits bifurcations. This result 
is further elaborated in Theorem 5.2, where we show that branches of formally stable 
relative equilibria satisfying a dimensionality condition generically have a continuous 
range of bifurcation points. This property is specifically related to the existence of 
continuous stabilizers and shows that these bifurcations are purely singular phenomena, 
since it is a well known result, originally due to Arnold [1], that branches of nondegenerate 
regular relative equilibria cannot bifurcate. In Section 6 the results of this article are 
illustrated in two study cases, the sleeping Lagrange top and the system of two point 
vortices on the sphere. 

We use a mixture of notation for derivatives: dh denotes the differential of the function 
h, and if /i : iV x X —>• K is defined on a product of spaces, Djyli denotes the partial 
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derivative(s) of h with respect to the A^-variables. Occasionally, we write d^/i to denote 
the Hessian matrix at the point z, but omit reference to the point z if it is clear. 

2. Hamiltonian relative equilibria 

In this section we provide necessary and sufficient conditions for the existence of 
Hamiltonian relative equilibria. The study will be local, in a neighbourhood of a given 
point in a symplectic phase space of a symmetric Hamiltonian system. To this end, we 
will use the framework provided by the bundle equations obtained in [23] and generalized 
in [24] (see also Chapter 7 in [16]). In our approach we will additionally incorporate in 
a critical way the freedom in several choices that exists if the symmetry group action 
exhibits continuous isotropy groups. 

To hx notation, we recall the standard definitions. A proper C?-Hamiltonian system is 
a quintuple {V, uj, G, J, h), where {V, uj) is a symplectic manifold, G is a Lie group acting 
properly and in a Hamiltonian fashion on V. We write the action by concatenation: 
(g, z) ^ g ■ z. The resulting momentum map 3 '.V ^ Q* satisfies 

= d(J(-),0. V^Gg. 

Here ^•p(z) = ^ _ pexp(t^) • z is the fundamental vector field associated to the element 

^ G g (the Lie algebra of G) evaluated at z £V. We assume J to be G-equivariant with 
respect to the coadjoint action on g*: 

J(g-z)=Ad;_i(J(z)), WgGG^zGV, 

where Ad* : G x g* — >• g* is the coadjoint representation of G. See also Remark 2.4 
for the adaptation to coadjoint actions modified by a cocycle. Finally, h G G°°('P) is a 
G-invariant Hamiltonian function on V. The symplectic manifold V is the phase space 
for the Hamiltonian dynamical system given by the flow of the vector field Xh defined 
by Hamilton’s equations 

LXh^ = 

Since h and u are G-invariant, the Hamiltonian vector field Xh is G-equivariant and 
so is its flow, therefore sending group orbits to group orbits for all times. A well known 
theorem by Noether states that J is constant along integral curves of Xh- Our main 
object of study is defined next. 

Definition 2.1. A relative equilibrium is a point in V whose integral curve lies in a 
group orbit. 

The following characterizations of relative equilibria are standard and can be found, 
for instance in [9]. 

Proposition 2.2. Let (V,u},G,3,h) be a G-Hamiltonian system and z G P with mo¬ 
mentum J(z) = g. The following are equivalent. 

(i) z is a relative equilibrium. 

(ii) The group orbit G ■ z consists of relative equilibria. That is, if z is a relative 
equilibrium, then so is g ■ z for any g in G. 

(hi) There is an element ^ G such that Xh{z) = ^'p(z), where Q^. < Q is the Lie 
algebra of the stabilizer of g defined by G^ = {g G G : Ad*-i/j = g}. 
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(iv) There is an element ^ such that z(t) = • z, where z{t) is the integral 

curve of Xh with z(0) = z. 

(v) There is an element C S such that z is a critical point of the augmented 
Hamiltonian 

h^ := /i-J«, 

with J^(z) := (J(z),^). 

Here and throughout, we use < to denote the relation of being a Lie subalgebra or a 
closed subgroup (according to context), and <1 to denote the relation of being a normal 
subgroup, or an ideal in a Lie algebra. We denote the normalizer of a subgroup H in G 
by Ng{H). 

Remark 2.3. The element ^ associated to z is called a velocity for the relative equi¬ 
librium. Notice that f is not uniquely defined if z has continuous isotropy, and this is 
a key observation that will follow a long way. The isotropy group, or stabilizer, of z, 
Gz = {g & G : <7 ■ z = z} is a compact Lie subgroup of G by the properness hypothesis 
on the action. If Gz has positive dimension it has a non-trivial Lie algebra Qz whose 
elements ly satisfy by definition rjviz) = 0. 

Notice also that the equivariance property of J implies in particular that g^ < g^. 
Therefore, by (iii) in Proposition 2.2 if ^ is a velocity for z, so is ^ -I- ry for any t] G Qz. 
The converse is also true: any two admissible velocities for z must differ by an element 
of g^. Also note that it follows again from (iii) and (iv) in the same proposition that if 
^ is a velocity for z, then Adg^ is a velocity for the relative equilibrium g ■ z. 

Remark 2.4. It was stated earlier that we are assuming the momentum map to be 
equivariant with respect to the coadjoint action g*. If G is compact, this can always 
be arranged by an averaging argument [10, 16]. However, for more general groups, J 
may not be able to be chosen to be equivariant in this sense. In such cases, there 
is a cocycle 0 : G ^ g* (the ‘Souriau cocycle’) for which the resulting modified action 
Coadg/i := Ad*g-ifi + 9(g) renders the momentum map equivariant: 3(g-z) = CoadgJ(z). 
Moreover, as shown in [25] there is an analogous local model for the symplectic action 
of G and similar bundle equations. It follows that all the results of this paper stated for 
the coadjoint action also hold in this more general case of a modified coadjoint action, 
provided the momentum isotropy group G^ is understood to be relative to this modified 
action. 

We state a well-known algebraic lemma. 

Lemma 2.5. Let H be a compact subgroup of the Lie group G, and let q — (B m be an 
Ad(H)-invariant decomposition (as vector spaces). 

(i) The subspace of fixed points g^ is a Lie subalgebra of g, and the Lie bracket on 

descends to one on m^. 

(ii) The Lie algebra of Nq(H) satisfies £,it(NG(H)) = f) © = fy + g^. Conse- 

guently, if L = Nq(H)/H, and I is its Lie algebra, then I ~ as Lie algebras. 

(iii) If H is a normal subgroup of G then = m. 

In fact = m is equivalent to the (weaker) statement that Ng(H) contains the 
connected component Go of G. 
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Proof. The decomposition follows by using an Ad(i7)-invariant inner product, which 
exists because H is compact, and defining m = 

(i) It is easy to check that is a Lie subalgebra of g. Now, t)'^ is in the centre of 
, and it follows that the Lie bracket descends to ~ 

(ii) For any f G g write ^ = + ^ f)©na. Now g G Na{H) means gHg~^ = H, and 

putting g = exp(t^) and differentiating at t = 0 shows ^ G n := 2,\z{Ng{H)) if and only 
if ^ satisfies the (linear) condition 

Ad^i^ — ^ G [}, V/i G H. 

Writing ^ = 17 + '?'*“ G l)©tTi (an Ad(iJ)-invariant decomposition), this condition becomes 

= 0 , 

which is equivalent to G m^, as required. 

(iii) This follows immediately from the fact that n = f) + rri'^. □ 

2.1. The bundle equations with isotropy. 

Since we are interested in the local properties of a Hamiltonian flow in a neighbourhood 
of a relative equilibrium, we will substitute the phase space V by the symplectic tubular 
neighbourhood given by the Marle-Guillemin-Sternberg (MGS) model. Originally due to 
Marie [8] and Guillemin and Sternberg [3], it is now standard and details can be found 
for instance in [16]. We briefly recall its construction. 

Let 2 ; G 7^ be an arbitrary point with momentum J{z) = g,. Let gz and g^^ be the 
Lie algebras of the stabilizers of z and g as before. We will assume through this paper 
that /j is a split element of g*, meaning that there is a G^-coadjoint invariant splitting 
g = g/i © q. Examples of split momentum values are the cases of compact or Abelian G^. 
For an example of non-split momentum see [24]. 

Since by the equivariance of J we have Gz < G^ (we use < to mean ‘is a closed 
subgroup of’), and since Gz is compact, we can find a Gz-invariant splitting 

g^=m©gz, (2.1) 

with associated dual invariant splitting g* = m* ® g*. Here g* and m* are identified with 
the annihilators of m and gz in g* respectively. Finally, combining these splittings of the 
Lie algebras, we have a Gz-invariant splitting g = m® gz ffi q. For each of these splittings 
we will denote by P[ the projection onto the factor 1. 

Next, choose N = A^z to be a complement to g^ • z in kerTzJ, Gz-invariant with 
respect to the induced linear representation of Gz on TzV. It follows from the general 
properties of Witt-Artin decompositions for Hamiltonian actions that such a choice is 
always possible. Moreover, {N, fl) is a symplectic vector space, where H is the restriction 
to N of w(z), and the linear action of Gz on N is Hamiltonian, with associated equivariant 
momentum map Jat : At —>• g* defined by 

{3N{v),rf) ^]^n[g-v,v), ^/r] £ Qz, v £ N, (2.2) 

where rj ■ v = ^ _ ge*** ■ v is the infinitesimal generator at v corresponding to 77 G gz- 

The space N is usually called the symplectic normal space (or symplectic slice) at z for 
the G-action on V. 
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The product space G x m* x iV supports free actions of both G and Gz , given by 

9'-{9,P,v) = i9'9,P,'>j) 9 ' (2.3) 

h-{9,P,v) = {gh~^,Adl-ip,h-v) h G Gz (2.4) 

for all g € G, p G m* and v G N. It is clear that these actions are free and commute. 
Consider the principal bundle associated to the G^-action 

7r:Gxm*xA^ -G- G Xq^ (m* x N) 

{g,p,v) HA [g,p,v], 

where [g, p, u] is the equivalence class of points consisting of the orbit Gz ■ (< 7 , p, v) C 
G X m* X iV. The G-action (2.3) on G x m* x iV descends to a smooth action on 
G Xa^ (nr* x N) given by 

9 ' ■ [ 9 , P, "y] = [ 9 '9, P, v] g' G G. (2.5) 

It follows from the MGS construction that there are open neighbourhoods of the origin 
in m* and N for which it is possible to define a local symplectic form on G xg, (nr* x iV), 
as well as a local G-equivariant symplectomorphism 

ip-.Gxa, (nr* xN)^V 

onto an open G-invariant neighbourhood of the orbit G • z, and (p satisfies (p([e, 0, 0]) = z. 
We will denote by Y the domain of this diffeomorphism in G Xg^ (nr* x N). The concrete 
expression of the mentioned symplectic form wy can be given as follows: Every tangent 
vector to Y can be written as ^\^p^i,{[g, p,v]) G p with 

1Kp,v{[9. P, I'D = T(g,p,v)T^{9 ■ A, p, h), 

where A G g, p G nr* and v G N. Here, g-X is the concatenation notation for Tf.LgX G TgG, 
where L : G x G —>■ G is the left translation on G, given by Lg/g = g'g. We then have 


wy (7 Ai ,pi ([p, Pw]), 7A2.P2 A2 ([ 5 , P, ■y])) 

= (p2 + T'„JAr(b2), Ai) — (pi + T„JAr(i;i), A 2 ) (2-6) 

+ (p + P + Jiv(ri)) [Ai, A 2 ]) + ^(hi, V 2 )- 

In addition, the induced G-action on (F, wy) given by (2.5) is Hamiltonian with a 
locally defined equivariant momentum map Jy=Jo(p:F ^g*. The explicit expression 
for Jy is very simple and provides in this local model a normal form for J. It is given by 

Jv([5,pw]) = Ad*-i(p-tp-f JAr(T)). (2.7) 

The idea exploited in [23] and [24] consists in taking the pullback by tp of the Hamil¬ 
tonian vector field Xh to F and then lifting it to G x m* x IV. Then the differential 
equations for the flow of a choice of lifted vector field are obtained, providing a general 
framework to study the local dynamics of Xh near G ■ z. 

In order to state these equations, note hrst that ho ip is a G-invariant function on F, 
so it lifts to a G X G^-invariant function h o tp o tt on G x m* x N. From equations (2.3) 
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and (2.4) it follows that 

{ho(po7T){g,p,v) = {(poTr)*h{g ■ {e,p,v)) = {(f o TT)*h{e, p,v) 

{hoTTO(f){l- {g, p, v)) = {(fio TT)*h{gr'^, Ad*-ip, I ■ v) = {(fi o TT)*h{e, Ad*-ip, I ■ v), 

for all g G G, I € Gz, p G vn* and v G N. Therefore we can identify h o ip o tt with a 
Gz-invariant function on m* x that we will denote by 

h e G°°(m* X (2.8) 

Since tt is a locally trivial fibration, a 7 r-projectable local vector field X G X(Gxm* xX) 
can be expressed as 

X(g,p,v) = {g ■ {Xg{g,p,v)),Xm-{g,p,v),XN{g,p,v)), 

where X^{g,p,v) G g, X,„t(g, p,v) G m*, and Xiq{g, p,v) G N. The bundle equations on 
G X m* X define a vector held on G x m* x A^ whose components are given by: 

Xg{g,p,v) = Dm*h{p,v) 

Xm-ig,P,v) = Pm- +JAr(t>))) 

XNig,P,v) = n'^{DNh{p,v)). 

Here Dm*h and D]s[h denote the partial derivatives of A S G°“(m* x N)'^‘ with respect 
to m* and N respectively, and H** : N* -G N is the linear Gz-equivariant isomorphism 
induced from fl. As shown in [23, 24, 16], these are the differential equations for the 
flow of the unique 7 r-projectable local vector held X on G x m* x A^ which is a lift of 
p*Xh G X(G xq ^ (m* X N)) (i.e. such that X and p*Xh are 7 r-related) and satishes the 
additional condition 

PgJXg)=0. (2.9) 

Notice that since Dm»h{p, v) G m, the above equations imply that Pq(Xg) = 0, which is 
a consequence of the condition on p to be split, and it is not related to the condition 
(2.9). Notice also that since is a G-equivariant symplectomorphism, p*Xh is actually 
Xhoip, the Hamiltonian vector held corresponding to the local G-Hamiltonian system 
{Y,ujy, G, Jy, h o p), which is a local model for {V,uj, G, J, h). 

However, in order to obtain all the projectable vector helds on G x m* x X that are 
TT-related to p*Xh one needs to include all vertical vector helds tangent to the 7r-hbres. 
These are, in view of (2.4), of the form 

X{g,p,v) = ( 5 -r 7 ,ad*p, -ry-u), 

where 77 : G x m* x X ^ gz is an arbitrary smooth map. In the following we will hx 77 to 
be an arbitrary constant in gz) since the vector helds obtained in this way generate the 
module of vertical vector helds. 

Using the facts that ad*p G m* (since m* is Gz-invariant), ad*(JAr(u)) G gj (since 
Jw(?^) G g*), and r] ■ V = Ht*(d(Jjv(-)i ^)) (since the Gz-action on N is Hamiltonian), we 
can write the equations of the how of the most general local vector held on G x m* x X 
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that projects to if*Xh as 

g = 9-iDm»Hp,v)+v) ( 2 . 10 ) 

p = Pm* + JAr(i;))^ (2-11) 

V = n\DNh{p,v)-d3l{v)). ( 2 . 12 ) 

where ?? S 0 z is arbitrary and the function J)(f is defined by 3]^{v) := (JAr(w),? 7 ) for all 
p G Qz and V G N. 

We also recall for future reference the following identity, which is satisfied for integral 
curves of the above flow (see equation (7.7.9) in [16]). 

Pgj (p + JAr(v))^ = TyJ]y{v). (2.13) 


2.2. Characterization of relative equilibria. 

We can now reduce the problem of finding relative equilibria of the G-Hamiltonian system 
{'P,uj,G,3,h) near z gV to find solutions of (2.10),(2.11),(2.12) that project to relative 
equilibria of Xhoip near [e, 0,0] in Y. We then have the following characterization. 


Proposition 2.6. A point [p,p, u] near z = [e,0,0] G Y is a relative equilibrium for 
Xhoip if and only if there exists an element p G Qz such that 


or equivalently 


adD„.^(p.„)+r,(P + Jw(^)) = 0. 

DNhnip,v) = 0, 


Pm* (^^^D^,h{p,v)+7J^P + ^n{v))^ 
DNhr,{p,v) 


0, and 

0 , 


(2.14) 

(2.15) 


(2.16) 

(2.17) 


where 

hr,{p,v) := h{p,v) - 

The element of q given by Adg{Dm*h{p,v) + p) is a velocity for the relative equilibrium 
[5,P, I']- 

Proof. Let [g,p,v] be in the domain Y of the diffeomorphism (p. First, by Remark 2.3, 
[g, p, u] is a relative equilibrium if and only if [e, p, v] also is, so we can restrict the 
study to points of this form. Let [e,p, u](t) be the integral curve of <p*Xh with initial 
condition [e, p, u](0) = [e, p, w]. Therefore, we can find representative curves g{t), p{t), v{t) 
with initial conditions p( 0 ) = e,p( 0 ) = p, u( 0 ) = v, respectively, such that [e,p, u](t) = 
[g{t), p{t),v{t)]. According to Proposition 2.2, [e,p, u] is a relative equilibrium for the 
local G-Hamiltonian system on G (>n* x N) corresponding to the Hamiltonian vector 
field ip*Xh = Xhoip, with velocity ^ G g if and only if 

t^QlGP,v]{t) =^v([e,p,u]). 


d 

dt 
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This is equivalent, in view of (2.3) to 


5(0) 

= £, + p' 

m 

= ad;,p 

u(0) 

= —p' • V 


for some rj' € gz- Since in the equations (2.10),(2.11),(2.12), rj is arbitrary, we can absorb 
r]' into rj so that the relative equilibrium conditions are 

^ = Dm*h{p,v) +r] 

0 = Pm* + JAr(i;))^ (2-18) 

0 = DNh^{p,v) 

Note that using (2.10),(2.11),(2.12) with this choice of isotropy rj, we have ^(O) = 0. 
Therefore, from (2.13) it follows that Pg* ^ad^^^ + JAr{i;))^ = 0. This, to¬ 
gether with (2.18) is equivalent to ^^_|_^(p-I-JAr(u)) = 0, since Dm*h{p,v) +r] G 

g^, p -I- Jjv(i^) S 0^ and g* = g* © m*. Therefore we have obtained that [e, p, v] is a rela¬ 
tive equilibrium ox\G (na* x N) near [e, 0,0] with velocity ^ if and only if there exists 
V G gz such that ^ = Dm*h{p,v) + p and Djs[hzf{p,v) = 0. And by Remark 2.3, this is 
equivalent to [g, p, i;] being a relative equilibrium with velocity Adg{Dm*h(p, v) + p). □ 

Remark 2.7. Note that (2.11) implies that the m*-components of the vector fields on 
G X m* X A' that are 7r-related to Xto^p are always tangent to the coadjoint orbit of 
Gg, which contains p + JAr(r')- This equation can be studied by topological methods, as 
in [10]. It can also be simplified by imposing conditions on G^. For instance, if G^ is 
Abelian, (2.14) is automatically satisfied. Also, the study of some dynamical properties 
can be simplified if these coadjoint orbits are bounded, as we shall see in Theorem 3.6. 
On the other hand, equation (2.15) is a critical point equation for a family of functions 
on N parametrized by (ly, p) G g^ x m*. Therefore, (2.15) is naturally suited to be studied 
by methods based on the implicit function theorem or singularity theory. Note also that 
in general hg G G°°(m* x N) is no longer G 2 -invariant, due to the presence of p, but 
only (G 2 ),;-invariant. 

2.3. Two lemmas. 

We will now focus on the solutions of equation (2.15) under the initial assumption that 
0 G A is a critical point of hr,{0, ■) for some p G gz- Notice that the dependence of on 
p is only through the term Jjf and this is a quadratic polynomial on A. It follows that 
if 0 is a critical point of hg{p, •) for some p G gz then 0 is also a critical point of h-q'ip, •) 
for any other element in p' G gz, since 0 is always a critical point of The problem 
consists then in finding triples {p, p,v) G m* x gz x N satisfying 

DNhg{p,v) = 0 , 


under the initial assumption 
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In order to attack this problem, we state two main lemmas, based respectively on non¬ 
degeneracy and degeneracy properties of suitable restrictions of a function on gz x m* x N. 
These two lemmas will be used in most of the results of this paper. 

Recall that if G acts on a set X, and K is a subgroup of G then one writes 

= {x £ X : k ■ X = X for all k £ K} 

for the fixed point set. Moreover, the action of G on X restricts to an action of the 
normalizer Na{K) on X^ , and hence defines an action on X^ of Na{K) / K . Before 
stating the lemmas, we recall the Principle of Symmetric Criticality, due to Palais [17], 
which we will have recourse to a number of times. 

Principle of Symmetric Criticality Suppose a Lie group G acts properly on a man¬ 
ifold M, and let f : M ^ M. be a smooth G-invariant function. Let x £ . Then 

df{x) £ (TfM)^, which as H is compact can be identified with In particular, 

it follows that if x is a critical point of f\j^H then it is a critical point of f. 

We now have the first of the technical lemmas. 


Lemma 2.8 (non-degeneracies). Suppose that H < Gz and let f £ G°°{QzXm* x N)^ 
be a smooth H-invariant function. Suppose there is a subgroup K < H, satisfying 

Dn f {0,0,0) = 0, and 
D%f{0,0,0)\ is non-degenerate. 

Then, there is a unique local (defined in a neighbourhood of (0,0)j smooth Nh{K)- 
equivariant map v : m*^ x satisfying u(0,0) = 0 and 


DNf{ri',p,v{p,r]')) = 0 


for every {p,ri') £ m*^ x near (0,0). 

The stabilizer of m = {p,v{p,ri')) £ m*^ x is {Gz)m 
satisfies 

K < {Gz)m < {Gz)p- 


{Gz)p n {Gz)y(^p^rj'^ and 


Proof. Notice that m*^, and support linear actions of Nh{K), which is the 
maximal subgroup of H that leaves invariant these subspaces. By the pull-back property 
of Hessians, it follows that D%f{0,0,0)\^K = (/|jyk)( 0, 0, 0). If this bilinear form 

is non-degenerate, it follows from the Implicit Function Theorem that there is a unique 
map v : m*^ x —>• defined in a neighbourhood of (0, 0) such that 

DNK{f\i^K){{v' , p,v{p,r]')) = 0 

for every (p, rj') in m*^ x g^, and these are all the points in m*'^ x gf^ x satisfying the 
equation (/|JY«^)(^^ PtO) = 0 near (0, 0, 0). By the invariance properties of f\jyK, it 
follows that the map v is X//(i4r)-equivariant. 

Since f{r]',p, ■) £ G°°{N) is in particular X-invariant for any {p,r]') £ m*^ x g^, it 
follows from the Principle of Symmetric Criticality, and v{p, rj') £ , that 


DNf{v',P,v{p,rj')) = 0. 


The property about the stabilizer of {p,v{p,ij')) is obvious. 


□ 
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The second technical lemma involves degenerate critical points, and will be applied 
to finding bifurcating branches of relative equilibria. It is usual to apply methods of 
singularity theory (e.g., [22, 11]) to study the nature of critical points when a function 
with a degenerate critical point is deformed. However, to avoid assumptions on higher 
order terms, and discussions of finite determinacy and versal unfoldings, we consider here 
cases that can be reduced to one variable, where one can give a general existence result. 

An action of a compact Lie group AT on a manifold M is said to be of cohomogeneity one 
if the orbit space M/K is 1-dimensional. Note that any representation of cohomogeneity 
one is necessarily irreducible. Some simple examples are the representations of the trivial 
group on K, SO(n) on R” and SU(n) on ~ C”. 

Let (t(u) be a continuous real-valued function on a topological space. We say cr crosses 
0 at uq if a{uo) = 0 and in any neighbourhood of ug, there exist ui,U 2 such that a{ui) > 0 
and a{u 2 ) < 0. 

Lemma 2.9 (degeneracies). Let W and N be representations of a compact Lie group 
K, and let f S C°°{W x N)^ be a smooth K-invariant function. Let L < K, and A a 
path connected open neighbourhood of the origin in . Suppose that 

D]sif{X^ 0) = 0 for all A G A. 

Let Af := kerZl^/(0, 0) H . Suppose Af ^ 0 and that the following two conditions are 
fulfilled: 

(i) the representation of Nx{L) on Af is of cohomogeneity one (and hence irre¬ 
ducible), and 

(ii) suppose that cr(A) crosses 0 at X = 0, where (t(A) is the eigenvalue of D'j^ffiX, 0) 
(thus cr(0) = OJ. 

Let = Af (B S be an Nx{L)-invariant decomposition. Then for each sufficiently small 

V € Af there is a X = Xy G A and s = Sy G S (not necessarily unique) such that 
Lt N f (.Xy , V, Sy^ 0 . 

The stabilizer of the critical point m = (Xy,v,Sy) is Km = K\^ n Ky which satisfies 
L<Km<NK{L). 

Proof. Notice that by hypothesis, for any fixed A G A the function /(A,-) G C^{N) 
is L-invariant. Therefore, according to the Principle of Symmetric Criticality, a point 

V G is a critical point of /(A, ■) if and only if the restriction of /(A, •) to N^, denoted 
by /^(A, •), has a critical point at v. This is equivalent to iAjv/(A, w) = 0. Notice 
also that the pull-back property for Hessians implies 

= D%Lf^{X,v). 

Consequently, if (i) and (ii) are satisfied, then the (unique) eigenvalue ct(A) ofZ?^j:,/^(A, 0)| 
is the only eigenvalue of Il^i/^(A,0) that crosses 0 at A = 0, the others being bounded 
away from zero. We now apply the Splitting Lemma from singularity theory, see for 
example [22, 11]. 

Let = Af (BS be the given splitting. Then Zl^/(0, 0)|g is non-degenerate. Writing 

V = {y,s) with y G Af, s G S, there is an A^/f (L)-equivariant change of coordinates of 
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the form (A, y, s) i—?• (A, y, s(A, y, s)), such that 

/(A, y, s(A, y, s)) = Q(s) + g{X, y). 

Here Q is the (L)-invariarLt non-degenerate quadratic form \D^f\g, and g is an 
(L)-invariant function for which D_\fg{X,0) — 0 and D'j^g{0,0) = 0. It follows that 
/(A, •) has a critical point at u = (y, s) G if and only if s(A, y, s) = 0 and y(A, •) has 
a critical point at y. 

Since we are only considering critical points with respect to N^, we can replace 
/(A,y, s) by /(A,y, s) — /(A, 0,0) and thus, without loss of generality, suppose that 
/(A, 0, 0) = 0. And hence y(A, 0) = 0. 

For the remainder, we first treat the case dim A/” = 1, and then extend to more general 
representations of cohomogeneity 1. With y G K and g(A, 0) = 0 and Dj>j-g{X,0) = 0 we 
can write g{X,y) = y^yi(A,y). It follows from (ii) that a{X) = 2yi(A,0). Therefore, by 
Taylor’s theorem, 

5 (A,y) = iy^cr(A) -Hy^5(A,y) 

for some smooth function b. Differentiating with respect to y we have 

gW y) = + y‘^c{x, y), 

for some new function c. Now use the fact that a crosses 0: let Ai,A 2 G A be such 
that cr(Ai) > 0 and a{X 2 ) < 0. There is then a neighbourhood V of 0 in AT such that 
(t(Ai) -I- yc(Ai,y) > 0 and a{X 2 ) + yc(A 2 ,y) < 0 for all y G H. Now fix y G H. By the 
intermediate value theorem there is a A G A (by taking any path in A joining Ai to A 2 ) 
with y'(A,y) = 0 as required. 

Now consider the case where dimA/” > 1, and let z{y) be the generator of the ring of 
invariants for the A^/f (L)-action on A/”. Then we can write g = y(A, z), and the argument 
above can be repeated using z in place of y. 

For the stabilizer, let m = (A, y, s) be a critical point. Then s = s(A, y, 0) (i.e., s = 0) 
so the result follows. □ 


3. Stability 

In this paper we use equations (2.10), (2.11), (2.12) together with (2.14) and (2.15) 
to obtain results on the nonlinear stability, persistence and bifurcations of Hamiltonian 
relative equilibria. To do so, we take advantage of the available freedom in the isotropy, 
given by the indeterminacy in the Lie algebra element y G 02 in the bundle equations 
with isotropy, and clarify the appearance in the literature of “orthogonal velocities” 
and optimal stability conditions. In this section we will be concerned only about the 
stability problem of relative equilibria. We start by introducing some notions needed 
in the following. We will always assume that we are given a G-Hamiltonian system 
(P, w, G, J, h), a point z with split momentum J(z) = y, and that we have chosen 
Gz-invariant splittings = 02 © in and ker r 2 J = 0 ^ ■ z © A/, and a G^-invariant splitting 
0 = 0/i©q- The results in this section extend the main result of [12] to allow non-compact 
G^, although we still require 0 to admit a G^-invariant inner product. 

Definition 3.1. Let z G P be a relative equilibrium and ^ G 0 ;_j a velocity. Then the 
projection = IPm(C) G m is called the orthogonal velocity of z. 
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Note that the orthogonal velocity of a relative equilibrium is unique once a splitting 
(2.1) has been chosen, since all velocities of z differ by elements of Qz- On the other 
hand, in general there is no canonical orthogonal velocity for z, unless the splitting (2.1) 
is unique (see Remark 3.7 for a counterexample). 

Definition 3.2. Let z & V he a. relative equilibrium with momentum J(z) = fji. We say 
z is nonlinearly stable (or stable modulo G^) if for every G^-invariant neighbourhood U 
of z there is an open neighbourhood O of z contained in U such that the integral curve 
of Xh through any point in O is contained in U for all time. 

This definition of stability for symmetric Hamiltonian systems was introduced in [18] 
where it is shown to be a natural one in this setting, since it has been noticed that the 
existence of a conserved momentum map can allow a drift of the Hamiltonian dynamics 
along the G^-orbits which makes the obvious choice of orbital stability too restrictive 
in the symmetric Hamiltonian scenario. We will see below that stability modulo G/z is 
related to the following notion of formal stability. 

Definition 3.3. Let z he a relative equilibrium with momentum J(z) = We say 
that z is formally stable if it admits a velocity C & such that is definite. 

The lemma below will be of great importance in the exchange of information between 
a G-Hamiltonian system {'P,uj,G,3,h) and its local model (T, wv, G, Jy, h o ip). In 
particular, it shows that the concept of formal stability is independent of the choice of 
symplectic normal space N. First we fix some notation. 

Definition 3.4. We say that an injective map z : W ^ V, where W is a neighbourhood 
of 0 in a vector space, is a parametrized branch of relative equilibria if for every w G W, 
the point z(w) is a relative equilibrium, and all the relative equilibria of the branch 
are inequivalent in the sense that they all belong to different G-orbits. Let z(0) = z 
with J(z) = p. The points of the branch z(w) are said to be of the same symplectic 
type (or symplectic orbit type) if for every w G W, there is some g € G such that 
Gj(^z{w)) ~ 9 ^and Gzf^w) ~ 9^z9 

Notice that by construction of the Marle-Guillemin-Sternberg model, if z{w) is a 
family of constant symplectic type then the symplectic normal spaces and Nz 

are isomorphic as G^-modules. Since group orbits of relative equilibria consist of relative 
equilibria, it is customary to refer, by extension, to G • z{w) as a parametrized branch of 
relative equilibria, although in this case we are not parametrizing points but G-orbits of 
relative equilibria. 

For the following result, recall that the function h is defined in (2.8). 

Lemma 3.5. 

(i) Let z = [e, p, 0] G Y, with p G m*. Then for every (p, p') G m* x Qz, 
DNhn'ip,-)=0 if and only if d[g_p^o] (^ o 


where 


ih o (p)^ := {hop) — Jy(-) = h^o p, 
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and ^ = Dm*h{p, 0) + 77 '. Furthermore, if 0 € N is a critical point of hr^'{p', ■), 
then 

Djfh^'ip,0) = 

where D^hjji{p,0) is the Hessian of the function hj^i{p, •) G C°°{N) at 0 € N 
and p 0 ] ° In restriction to N of the Hessian of {h o ip)^ at [e, p, 0]. 

(ii) Let z £ V be a critical point of h^ for some f £ Q- Suppose that is 

non-degenerate, degenerate, positive definite or negative definite. Then the cor¬ 
responding property is satisfied for any other choice of symplectic normal space 
N'. 

(iii) The same conclusions hold if in (ii) we replace N by and N' by {N')^ 
respectively, for any compact subgroup K <Gz. 

Proof. We will assume, without loss of generality, that Gjy ([e,p,o])) = Gp for each point 
of the family. In particular, from (2.7), this assumption implies that (G'p)p = Gp- Recall 
that from (v) in Proposition 2.2 the critical points of {h o ip)^ correspond precisely to 
relative equilibria of the G-Hamiltonian system {Y,ujY,G,JY,h o p). Now using (2.7) 
and 

{hop){[g,p,v\) = h{p,v) 

we have that 

(d(/i o v?) 4 ([g, p, u]), 7A,p.i,([5, p, u])) = 


It follows that if we make g = e and v 
93 )^([e, p, 0 ]) = 0 if and only if 

Dm-h{p,Q )= 0 
7?A,()i-(Jjv,P0jO))(p,O) = 0 

adj (p + p) = 0 

Therefore since Pgj(‘C) = p' the first condition is automatically satisfied and the second 
is equivalent to 

D]sihr,'{p,Q) = 0 . 

Notice that since ^ G 0 p and (Gp)p = Gp the third equation is automatically satisfied. 
This proves the first part of (i). 

For the second part of (i), notice now that if 0 G is a critical point of hrj'{p, •)> then 
D]Yh{p,0) = 0 and D3jf{0) = 0 for any p G m*. This follows from (2.2) since is a 


dt 


t = 0 
d 


{h{p -\-tp,v-\- tv) 


dt 


j _ p(Ad*ggtA)-i(p + p + tp + JAr('C + tv)), f) 

Dm* h{p, v) ■ p-\- DNh{p, v) ■ V 

+ (Ad*-i (adA(p + p + JAr(T))), Cl 
-{Adl-.{p + DM3N{v)-v),C 
= 0 in the previous expression we obtain d(/i o 
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quadratic polynomial on N. We can naturally identify N with the subspace of T'[e,p,o]^ 
given by { 7 o,o,{)([e, p, 0]) : v € N}. Then, we have 

° (7o,o>i ([e, p, 0]), 70.072 ([e, P, 0])) 

= Djfh{p,0){vi,V2) - T>^J^(0)(hi,i>2) 

= D%hr,>{p,0){vi,V2) 

with ^ = Dm* h{p, 0) + p'. 

For (ii), let N and N' be two symplectic normal spaces at z. That is, both are Gz- 
invariant complements to ■ z in ker T^J. Every element in N' can therefore be written 
as i; + A • z for some v € N and A S Since the Hamiltonian vector field is equivariant, 
the point p • z is a relative equilibrium with group velocity Adg^, and hence for all g G G, 
the corresponding differential vanishes: dhAdg^ig • z) = 0. Write g = exp(tp) for p G g^, 
and differentiate with respect to t at t = 0 to obtain 

d^/i5(p-z,-) -dJ[^_5] = 0, 

where the differentials are taken at z. It follows that for any v G ker dJ(z) and any p G g, 
we have 

d^/ij(p ■ z, v) = 0. 

It then follows that, given any A G g^ (so that A • z G kerdJ(p)) 
d^h^{v + A • z, i; + A • z) = d‘^h^{v, v), 

as required. Part (iii) is a straightforward consequence of the pull-back property of 
Hessians. □ 

We can now extend a standard stability result, originally stated in [5] and [15] for 
the case of relative equilibria with continuous isotropy, extending the work of [18] on the 
nonlinear stability of regular relative equilibria. This also extends the theorem in [12], 
where it was assumed that Gp is compact. 

Theorem 3.6. Let {'P,ui,G,3,h) be a G-Hamiltonian system and z G V a relative 
equilibrium with momentum J(z) = p. Suppose that g admits a Gg^-invariant inner 
product with respect to the adjoint representation. Suppose in addition that z admits a 
velocity ^ G g^ such that d^h^^^ is definite, where N is any Gz-invariant complement to 
Qg- z in kerT^J. Then z is nonlinearly stable, in the sense of Definition 3.2. 

Proof. Since (p : Y = G (nr* x N) V is a local G-equivariant symplectomorphism, 
it is clear that [e,0,0] is a relative equilibrium with momentum Jv([e,0, 0]) = p for 
the local G-Hamiltonian system {Y, cuy, G, Jy, h o p) that models the original system in 
a neighbourhood of the group orbit G ■ z. It is easy to verify that N, identified with 
{7o.o.i([e,0,0]) : V G N} G T[,,o.o](G Xg. (m* X N)) is a symplectic normal space at 
[e,0, 0]. Let C ^ g^ be a velocity for [e, 0,0] and let ^ -I- p G m 0 g^ according 

to the splitting (2.1). Then, it follows from (i) in Lemma 3.5 that Z3Ar/i,,(0, 0) = 0 and 
Dm*h{0, 0) = ^■*‘. Suppose now that d^^ ^ (h o is definite. Then, from (ii) also in 

Lemma 3.5 we have that Djghrj{0,0) is definite. The existence of a Gp-invariant inner 
product on g implies that p is split, so (2.10-2.12) apply. 
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By the smoothness of the dependence of h on m*, it follows from the Morse lemma 
that 0 G is a Lyapunov stable point for (2.12). Also, the condition on the existence 
of the inner product guarantee that the coadjoint orbits of are contained in compact 
hypersurfaces, so they are bounded. This, using (2.11) implies that 0 G m* is Lyapunov 
stable for (2.11). Since Dm*h{p,v) +?? G for all {p,v) &m* x N, it follows from (2.10) 
that g(t) G goGp, for any initial condition g(0) = go and all t. 

Let Og C G be an open neighbourhood of e. Any G^-invariant neighbourhood of 
[e, 0,0] in G Xg, (m* x N) is given hy U = {G^ ■ [g, p,v] : g e Ug, p & !/„•, v G Un} C 
G Xg^ (m* X N) where Ug C G, Gm* C m* and Un C N are neighbourhoods of e,0 and 
0 respectively. It follows from the above discussion, and from the Lyapunov stability of 
(0, 0) in m* X A^ that we can find open neighbourhoods Om* C Um* and On C Un of 
the origins such that p{t) G Um* and v{t) G Un for all t if p(0) G Om* and if z;(0) G Oat. 
Therefore, calling O = {[g,p,v] : g G Ug, p G Om*, v G On} C G Xg^ (m* x N), we 
have that the integral curves of Xho,p through points in O always lie inside U, proving 
the nonlinear stability of [e, 0,0], hence of 0 . It follows from Lemma 3.5(ii) that this is 
independent of the choice of symplectic normal space N. □ 

Remark 3.7. Results analogous to Theorem 3.6 have been obtained first in [18] for 
the case when G^ is discrete and in [5, 15] under similar hypotheses as here. In [5] the 
result is proved using the symplectic slicing technique, and the proof of [15] is based on 
an extension of the methods of [18] to the case of continuous isotropy. Theorem 3.6 is, 
however, more general, since it could guarantee nonlinear stability in some cases in which 
the application of [5, 15] is inconclusive. The main difference with our result is that in 
those two references, the bilinear form on N used to test stability is dz/ijr, where 
is the orthogonal velocity associated to some choice of Gz-invariant splitting (2.1). The 
strongest stability results follow then by testing over all possible orthogonal velocities 
corresponding to all possible such invariant splittings. In our approach, we fix from the 
beginning one splitting, which is necessary only for the construction of the local model 
of the G-Hamiltonian system, and we test over all possible velocities admissible for the 
relative equilibrium under study. This is the same as testing stability for + 77 for all 
possible 77 G 0z, where is the orthogonal velocity corresponding to the fixed invariant 
splitting of Qfj,. We provide now an example where the methods of [5, 15] don’t predict 
stability, but Theorem 3.6 does. This is a different example to the one in [12], and here 
the relative equilibrium is not an equilibrium. 

Let z be a relative equilibrium with = SO(3) and Gz = S^. If we identify 
with and gz with span(e 3 ), the only Gz-invariant complement to gz is given by 
m = span(ei,e 2 ), which can be identified with m* using the standard inner product on 

Suppose that the symplectic normal space A^ at z is isomorphic to R^, where in a 
Darboux basis the symplectic matrix 12 takes the form 

/ 0 10 0 \ 



\ 0 0 - 10 / 
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Suppose that in the linear coordinates {xi,yi,X 2 ,y 2 } associated to this basis we have 

Kp,v) = ^{xl+yl-xl- yl) + /(p), 

where v = (xi,yi,X2,1/2), p G m* and / is an arbitrary S'^-invariant function on m*. 
Suppose also that in this basis for N, the symplectic action of Gz = is given by 




Re 0 
0 Rg 


with Re the standard rotation through 0 on The associated momentum map has the 
expression 

Jn(w) = ^{xl+yl+xl+ p|) . 

It follows from Lemma 3.5 that testing definiteness of with being the orthogo¬ 

nal velocity relative to the invariant splitting of into gz and m is equivalent to testing 
definiteness of 

i5po(0,0) =diag(l,l,-l,-l), 

and the test is inconclusive. Since there are no more invariant splittings available, the 
results of [5, 15] can’t predict the nonlinear stability of this relative equilibrium. In our 
approach, using Theorem 3.6 it follows from Lemma 3.5 that the relative equilibrium will 
be nonlinearly stable if there is some p G R ~ 02 for which Djfhrj{0,0) is definite. We 
have the general expression 

S(p, v) = {xl + yl) - {xl + yl) + f{p). 

Then we obtain 


D%hr,{0, 0 ) = diag(l - p, 1 - p, -(1 p), -(1 p)), 

which is definite for p G {— 00 ,—l)U(l,-|-oo), and therefore the relative equilibrium is 
nonlinearly stable. 

Finally, we remark that the first time that an approach similar to ours was used in 
the context of testing nonlinear stability of Hamiltonian relative equilibria is in [26]. 
In that reference, the authors only consider free Hamiltonian actions, which do not 
exhibit isotropy, however in that situation they investigate more general momentum 
values without requiring the existence of a C?^-invariant inner product on g. 

4. Persistence 

The problem of persistence of relative equilibria consists in providing conditions under 
which a given relative equilibrium z belongs to a continuous family of relative equilibria, 
called the persisting branch. Additional properties about this family, such as the stability 
and stabilizers of its elements, or its geometric features, are usually also of interest. 

Definition 4.1. A relative equilibrium z is said to persist, if for every G-invariant 
neighbourhood U containing z, the set U\G • z contains a relative equilibrium. 
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In several references, it is also required in the definition of persistence that the per¬ 
sisting set of relative equilibria is the G-saturation of a parametrized branch of relative 
equilibria of the same symplectic type, as in Definition 3.4. We will not require this. In 
this section we will examine four different scenarios where a given relative equilibrium 
can persist, together with additional information about the persisting branch. 

4.1. Persistence in the case of non-degeneracies. 

The following theorem gives sufficient conditions for finding persisting branches of rel¬ 
ative equilibria starting with a relative equilibrium satisfying a typical non-degeneracy 
hypothesis, as well as an estimate on the orbit types of the persisting relative equilibria. 
It extends Theorem 3.2 of [23] to non-compact groups, and uses a weaker hypothesis. 

First recall that the orbit type of a point z € M is the conjugacy class in G of its 
stabilizer Gz and it is denoted by (G^). In particular, two distinct relative equilibria 
have the same orbit type if and only if they have conjugate stabilizers. 

If M is a smooth manifold on which the Lie group G acts smoothly and properly, for 
each orbit type (H), the set 

’^(H) = {z G M : Gz is conjugate toH}, 

called the manifold of orbit type (H) in V, is a disjoint union of connected embedded 
submanifolds of V. Our main persistence result in the case of non-degeneracies is the 
following. 

Before stating the theorems, we discuss briefly an algebraic hypothesis we make below 
and in other theorems: 

Lemma 4.2. Let H be a compact subgroup of the Lie group G, with Lie algebras f) < 
respectively. Suppose one can write g = [} © m (as vector spaces), with m C 3 ( 0 ), the 
centre of q. Then 

(i) the decomposition is -invariant; 

(ii) The decomposition g = () © m zs on equality of Lie algebras, with m Abelian. 

(iii) [ 0 , 0 ] C f); 

(iv) any a,d{\))-invariant decomposition (in particular if it is Ad{H)-invariant) g = 
() © m' also satisfies m' c 3 ( 0 ); 

Furthermore, if G is compact then m C 3 ( 0 ) is equivalent to [ 0 ,g] C f). 

When this condition that m C 3 ( 0 ) is satisfied, we say f) is a co-central subalgebra of 
0 . It holds for any subalgebra if g is Abelian, and a non-Abelian example is provided by 
iL = SU(n) as a subgroup of G = U(n). 

Proof, (i) For the decomposition to be ad(f))-invariant simply means [f),m] C m, but 
since m C 3 ( 0 ), we have [(), m] = 0 c m. 

(ii) This is clear, since = [m,m] = 0. 

(iii) This follows immediately from (ii). 

(iv) Since g = ()©m = f)©m', the projection g —m with kernel () induces an isomorphism 
/ : m' ^ m. Moreover, since m' is ad([})-invariant, the projection commutes with ad(f)), 
so the two representations m and m' are equivalent. Thus [(),m'] = 0. Moreover, since 
m C 3 ( 0 ) it follows that [m,m'] = 0. Therefore [ 0 ,m'] = 0, whence m' C 3 ( 0 ), as required. 
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For the final part, suppose G is compact, and [g, 0 ] C f). It follows that f) < 0 (an 
ideal) since [g, ()] C [ 0 ,g] C h, and hence (by compactness of G), there is a decomposition 
0 = [} © m that is ad( 0 )-invariant. In particular, this implies [g,m] C m. But since 
[ 0 ,m] C 1) by hypothesis, we conclude [ 0 ,m] = 0, whence m C 3 ( 0 ). □ 

Note that since the adjoint action is by Lie algebra isomorphisms, fixed point sets in 
0 are necessarily subalgebras (see Lemma 2.5). In particular if iL < G^ then is a 

subalgebra of g^, and moreover (since K is then compact), (g*)^ is naturally the dual 
of 0 ^. The lemma above will be applied with g^ in place of () and g^ in place of g. 

Theorem 4.3. Let z G V be a relative equilibrium with momentum 3{z) = and let 
^ G Qfj. be a velocity for z. Write ^ = C‘'‘ + ? 7 Gm© 0 z according to the splitting (2.1). 
Suppose there exists a subgroup K < (G^),, such that 

(i) is a co-central subalgebra of , and 

(ii) dfh^\j^K is non-degenerate. 

Then, there is a smooth N(^Q^-j^(K)-equivariant map 

z : X 0 z^ —^ 

(P, v') ' — z{p, 77 ') 

defined in a neighbourhood of ( 0 , 0 ) such that z{0, 0 ) = z, and for each p' the map 
p i-A z{p, p') is an immersion, such that for each {p, p') in the domain, the point z{p, p') 
is a relative equilibrium with velocity of the form ^ -\- p' G and stabilizer Gz{p,r]') 
satisfying 

K < Gz{p,r]') < (Gz)p. 

Moreover the branch z{p,p') consists of every relative equilibrium near z with stabilizer 
containing K and velocity satisfying ^ G g^. 

This theorem if only non-trivial if m*^ 7 ^ 0. This is because, if m*^ = 0, then the 
theorem is satisfied by the constant map z{p') = z, since z is by assumption a relative 
equilibrium with velocity ^ p' for every p' G Qz- 

Proof. Write ^ = ^'*‘+77em©0z according to the splitting (2.1), thus p := Pg,,(^) G Qz- 
Using the local model (F, wy, G, Jv, h o ip) around the group orbit G • z, and putting 
f{p', p, v) = h^+r)' {p, v) the hypotheses imply that DNf{0, 0 ,0) = 0 and Z3^/(0, 0,0) 
is non-degenerate (see (2.8) for definition of h). Therefore, from Lemma 2.8 it follows that 
there exists a local smooth (iL)-equivariant map v : m*'^ x g^ satisfying 

Dn f{p', p,v{p,p')) = DNhr,+r]'{p,v{p,p')) = 0. Therefore (2.17) is satisfied for the 
isotropy Lie algebra element p + p' . In the local model, we therefore put z{p,p') = 
{p,v{p,p')), which for each p' is clearly an immersion. We will now show that (2.16) 
is also satisfied for the same choice of element. Using the Gz-invariant (and therefore 
if-invariant) splittings g^ = m©gz and g* = m*®g*, and noting that p-\-3 n{v{p,p')) G 
m* © gj = g* we can identify p + 3m{v{p,p')) with an element A G g*. Moreover, 
A e since p G and by the equivariance of Jat and the fact that v{p, p') G , 
3n{v{p,p')) G g*^. Therefore, (2.16) is equivalent to 

Pm* (ad* A) = 0 , 


(4.1) 
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with 7 = Dm* h{p, v{p, 'ri'))+'r]+r]'■ We claim that 7 S . It then follows from hypothesis 
(i) that the equation above is satisfied. To see this, note that (4.1) is equivalent to 

(ad* A, 13) =0, V/3 G m^. 

This in turn is equivalent to (A, [ 7 ,/?]) = 0 (for all f3 G m^). The hypothesis implies 
C 3 ( 0 ^), whence [P,j] = 0 for all [3 G 7 G showing (4.1) to be satisfied, 
independently of the values of 7 G g^ and A G g))^- 

To prove the claim that 7 G g^, recall that K < (Gz)ri, and therefore p + p' G It 
remains to show that Dm*h{p,v{p,p')) G m^. Now h is G^-invariant, and therefore for 
each V G , the function p 1 —!> h{p,v) is invariant. It follows (from the Principle of 
Symmetric Criticality) that for p G m*'^ we have Dm*h{p,v) G as required. 

Since we see that (2.16) and (2.17) are simultaneously satisfied, if we define z{p, p') G V 
by z{p,p') = (p{[e, p,v{p,p')]) then it follows by the equivariance of p, Proposition 2.6 
and the fact that tp is a local isomorphism of symmetric Hamiltonian systems that for 
each pair (p, p') G m*^ x g^, the point z{p, p') G P is a relative equilibrium with velocity 

^ + v' ■ 

The stabilizer of z(p, p') is Gz[p^r)')- Since p is G-equivariant this is equal to G[e p = 

{Gz)m where {Gz)m is the stabilizer of m = {p,v{p,p')) G m* x N with respect to the 
diagonal G^-action. According to Lemma 2.8, K < {Gz)m (so z G V^) and is contained 
in {Gz)p- It also follows from Lemma 2.8 that the pairs {p,v{p,p')) G m*'^ x give 
all the solutions to equation (2.15) with fixed ^ and p, from which the last statement 
follows. □ 

4.2. Persistence of branches of constant orbit type. 

As an application of Theorem 4.3, we now study in more detail under which conditions 
we can guarantee that a relative equilibrium belongs to a smooth branch of relative 
equilibria of the same orbit type (and possibly the same symplectic type), parametrized 
by momentum values. The following is an extension within our framework of the results 
of Patrick [19] in the free case and Lerman and Singer [5] in the singular case, each 
reference requiring slightly different hypotheses. It is also related to results of Patrick 
and Roberts [20], where G is assumed to be compact and acting freely, and they give, 
inter alia, necessary and sufficient conditions for the set of relative equilibria to form 
a symplectic submanifold. It would be interesting to adapt their approach to non-free 
actions. 

Theorem 4.4. Let z GV he a relative equilibrium with momentum p = J(z) for which 
is non-degenerate for some admissible velocity f G Qp. Suppose moreover that 
the Lie algebra I of L := Nc,,^{Gz)/Gz is Abelian, where Ng,,^{Gz) is the normalizer of 
Gz in Gp. 

(i) There is a G-invariant neighbourhood U of z in 7^(0,) such that the set Z of 
relative equilibria in U forms a smooth submanifold of dimension dimG/Gz + 
dim))!!*)*^*. For each p' G J(G) there is a relative equilibrium z' G Z with 
J(z') = p'. 

(ii) Suppose in addition that Gz <3 Gp, or more generally Zie{NG,^{Gz)) = Qp- Then 
the submanifold Z is a symplectic submanifold ofV. 
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(iii) IfGfi is compact and Abelian then the momentum values of the relative equilibria 
of Z have stabilizers eonjugate to G^,, and the family Z is of eonstant symplectic 
type. 

(iv) Finally, still assuming the connected component of G^ to be a torus, suppose z 
is formally stable. Then all the relative equilibria of this branch close enough to 
z are also formally stable and therefore nonlinearly stable. 

Proof, (i) In the local model Y = Gxq^ (m* xN) oiV near G-z, the point z corresponds 
to [e,0, 0]. According to the property 

G[g^p^v]= 9{i.Gz)pF{Gz)v)g 
the points of orbit type (G^) are those of the form 

yg &G, p v & . 

We wish to find solutions of (2.14, 2.15) in As usual, the velocity f can be 

decomposed as ^ = ^'^ + g & m (B Qz, see (2.1). In the hypothesis of the theorem, we 
have that DNhrj{0,0) = 0 and Zl^h,,(0,0 )|^g2 is non-degenerate. Using that Jy is a 
homogeneous quadratic polynomial on N and that ry-u = 0 for every g £ Qz and v £ , 

we have that DNh{0,0) = DNhg{0,0) = 0 and that i4^/i(0,0)= U^h,,(0,0 )|^g* is 
non-degenerate. 

By the Principle of Symmetric Criticality, and the pull-back property of Hessians, this 
means that the function h^^{p,-) := h{p, with p G has a non-degenerate 

critical point at 0 G for p — 0. Therefore, from the implicit function theorem there 
is a smooth function u : —)• such that for any g £ G and p G the points 

in the branch [g,p,v{p)] are all the solutions of (2.15), with g = 0, with orbit type (Gz) 
for p sufficiently close to 0. 

At points of this branch, condition (2.14) is equivalent to ad^ ,h{pv{p))P ~ again 
because JAr(^^) = 0 for v G 

To see that this equality is satisfied for all p in a neighbourhood of 0, notice that since 
p G and v{p) £ then by the Principle of Symmetric Criticality, Dm*h{p, v{p)) £ 
(as in the proof of Theorem 4.3). Since m®"', together with the bracket induced 
from g is isomorphic to the Lie algebra I, and is isomorphic to (m*^*)* = [*, the 

hypothesis on the Abelian nature of 1 implies (2.14). 

By (2.7), and since 3n{v) = 0 for u G the values of JL near p form a 

neighbourhood of p in the set 

{Ad*-i(p -I- p) : g £ G, p £ 

Let p' be one of such values, close to the coadjoint orbit through p. Then, there are 
g £ G and p G close to 0, such that p' = Ad*-i(p -I- p). It follows that the point 

[g,p,v{p)] is a relative equilibrium of orbit type (Gz) close to the orbit G ■ z. 

Furthermore, if p G (g^*^" C m* C (m*)'^"’, then p' has orbit type (G^) and the 
resulting relative equilibrium will have the same symplectic type as z. 

(ii) The symplectic form on Y is given by ujy in (2.6). Pulling this back to the sub¬ 
manifold Z of relative equilibria in (i) gives a closed form. To show it is non-degenerate 
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in a neighbourhood of [e,0, 0] it suffices to show it is non-degenerate at that point. Let 
71,2 be arbitrary tangent vectors at this point: 

71,2 = Tnig ■ Ai, 2 , Pi, 2 , Dv{p) ■ pi, 2 ) S 

Note that pi ^2 G (m*)*^*, but by Lemma 2.5(iii) the hypothesis Gz<Gp implies = m, 
and so 

TzZ ~ g • z X = g • z x m*. 

It follows from (2.6) that at [e,0,0], 

‘^y(7i)72) = (P2, Ai) - (pi, A 2 )-I-[Ai, A 2 ]) 

-f nN{Dv{p) ■ pi,Dv{p) ■ P 2 ), 

since = T„(p)JAr|^G* = 0. We can decompose an element A € g as 

A = A®^ -f A"" -h A®, 

with respect to the fixed splitting g = © m © q. Then 

‘^y(7i)72) = (P2 j A)") — (pi, A^) + (p, [A^, A|]) 

+ nN{Dv{p) ■ pi,Dv{p) ■ P 2 ). 

The first two terms are the standard symplectic form on m x m*, the last term can be 
seen as a ‘magnetic form’ on the same space, which does not alter the non-degeneracy, 
while the third term is the standard KKS symplectic form on g° = q*. It follows that 
UJY is indeed non-degenerate on TzZ ~ g • z x m* = (m x q) x m*. (See also Lemma 4.2 
of [5].) 

(hi) If Gp is compact then, using Palais’ tube theorem, one can choose a G^-invariant 
slice to the coadjoint orbit through p, which can be identified with g*. Then the stabilizer 
of Jy([e, p, v{p)]) = p. + p is {Gp)p, and since Gp is Abelian, {Gp)p = Gp. Thus all points 
of Z have the same symplectic type. 

(iv) Finally, suppose z is formally stable; that is, |is definite for some admissible 
velocity ^ G gp. Since by (iii) the points of Z are of the same symplectic type, they 
have isomorphic symplectic normal spaces N' ~ N. By continuity, the appropriate 
nearby hessian is also positive definite, and hence by Theorem 3.6 the relative 

equilibrium is nonlinearly stable. □ 

4.3. Persistence in the case of degeneracies. 

We now prove a persistence result specific for relative equilibria with continuous isotropy, 
which exploits the indeterminacy of the velocity in order to predict branches of persisting 
relative equilibria. 

Theorem 4.5. Let z be a relative equilibrium with momentum J(z:) = p and suppose 
that z admits the velocity ^ G gp with ^ = ^^ + r]Gm(Bgz. Let K = {Gz)ri and suppose 
that there exists a subgroup L < K for which M := kerd^hfl^i, is non-zero. Suppose 
moreover the following conditions are satisfied 

(i) g^ is co-central in gp, 

(ii) the action of Nk{L) on N is of cohomogeneity one, and 

(iii) the eigenvalue cr{g') of d^/i^+r,' for 77 ' G g^, crosses 0 at p' = 0. 
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Then, for each v G A/’\{0} close to 0 there are elements r]'„ G close to 0 and Zv GV 
close to z which are relative equilibria with velocity f + ijy and stabilizer Ky, which 
contains L. 

In Section 7, we give an elementary illustration of this theorem in the system of pairs 
of point vortices on the sphere. 

Proof. The proof is a direct application of Lemma 2.9, where we choose W = , and 

f{r]',v) = hfj+r]'{Q,v) (with p = 0). Notice that if z corresponds to [e,0, 0] in the local 
model, then Lemma 2.9, predicts that for each v G J\f there is an such that n is a 
critical point of h^+ri'„{0, •)• Therefore, (0,n) satisfies (2.17) for rj + rjy. By an argument 
similar to the one used in Theorem 4.3, the Principle of Symmetric Criticality guarantees 
that Dm*h{0,v) + ri + r]y G gji, and then condition (i) implies that (2.16) is also satisfied. 
Therefore, in the local model the persisting relative equilibria predicted by this theorem 
are group orbits of points of the form [e, 0, v] for each v G M. These are the points that 
correspond to Zy in P. 

The stabiliser is given by {Gz)v C (Gz)?;'- But p' G gf, so {Gz)n' = K, and v G 
whence L < {Gz)v D 

4.4. Persistence from formally stable relative equilibria. 

To end this section, we prove a more general result on momentum parametrized branches 
persisting from formally stable relative equilibria. It was originally obtained in [13] 
in a slightly different setup, but there it was required that be compact (see the 
remark below for more details). Here we present it as a product of the bundle equations 
formalism. 

Theorem 4.6. Let z be a formally stable relative equilibrium with p = 3{z). Suppose 
that 0 ^ has a G^-invariant inner product (for example, G^ is compact). Then there is a 
G-invariant neighbourhood U of z such that for every p' G J(U) near p there is a relative 
equilibrium z' near z with J(z') = p'. 

Proof. Let Uq be a neighbourhood of z for which the Marle-Guillemin-Sternberg normal 
form is valid. Recall that for a point [g,p,v] GY := G (m.* x N) we have 

^{[9,P,v]) = Ad*g-i{p + p + 3n{v)) 

Therefore J(T) = {Ad*-i(/r + p + e) : g G G, p G m*, e e im Jat}. Let lb C T be the 
subset corresponding to Uq; then J(To) C J(T). If G^ is compact then J(To) is an open 
subset of 3{Y), but it is not known in the more general setting here—see also Remark 
4.7 below. 

Step 1. Taking p' = p + p + e 

Suppose z' G Yq (or Uq) is a relative equilibrium with J(z') = Ad*-i {p + p + e). This is 
equivalent to g~^-z' being a relative equilibrium with momentum Ad*(Ad*-i (/x+p+e)) = 
p + p + e. Therefore in the statement of the theorem we can choose, without loss of 
generality, p' = p + p + e. 

Step 2. Using the bundle equations 

According to Proposition 2.6, a point [g,p,v] is a RE if and only if there exists p G gz 
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such that 

^‘^*D„,Hp,v)+^iP + ^Niv)) = 0 (4.2) 

DNhrj{p,v) = 0. (4.3) 

Equation (4.2) is an equation in g* and is equivalent to 

+ = 0 (4-4) 

where O = {Ad*-i (/9 + JAr(u)) : I £ G^} C g* is the G^-orbit oi p + Jiv(t) for the 
coadjoint action. 

Now we need to use the formal stability of z = [e, 0, 0]. By Lemma3.5(i), this is 
equivalent to there being a 7 £ 02 (which we fix throughout) such that is definite 

(as a quadratic form on N). Let us assume it is positive definite, the negative definite 
case being similar. 

Define / £ G“( 0 * x N) by 

f{a,v) := h^{a\^,v) - h^{0,0) 

and let fo be its restriction to G x TV. Define also ^ : g* x TV —>■ g* by 

(l){a,v) := 3 n { v ) - al , 

I 02 

and (j)o its restriction to G x N. 

Now (4.3) and (4.4) are together equivalent to the restriction of fo to (()q^( 0) having 
a critical point at the point (p © JAr(z;), u) £ m* © g* x TV = g* x N, where p £ 02 is the 
Lagrange multiplier. Indeed, the latter can be written 

DnU - if), v)) = 0 
4^0j(/- 

where a = p +JAr(u) (that is, (j){a,v) = 0). The first of these is equivalent to (4.3), while 
the second becomes (4.4), since (/ — {(j), rj)) is independent of the g* component of a. 

Step 3. Existence of critical points 

We claim that for sufficiently small G, fo is proper. To see this, we can apply the 
Splitting Lemma from singularity theory, see for example [22, 11]: since D^f^z) = 0 and 
D‘j^f{z) is non-degenerate, there is a neighbourhood U of (0,0) in g* x and change of 
coordinates on U of the form (a, v) 1 —>■ (a, E), where V = V(a^ u), and a smooth function 
p : g* ^ R such that 

/(a, v) = g{a) + /(0,0) + D%f{V, E), 

and similarly for fo- For simplicity, let us assume /(0,0) = 0. Now let / = [a, 6 ] C R be 
a compact interval. Then, since G is compact, 

foia,v)Gl ^ 11^/(14,14) £ [a-maxp(G), 6 -minp(G)]. 

Since 11^/(14,14) is positive definite, this implies 14 in contained in a compact set TV/. 
Then 

fo\l)cOxNi 
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which is compact, and hence indeed fo is proper. It is also clear that fo is bounded 
below, by vaingiO). 

To finish we use a well-known variational argument: any continuous function f : X ^ 
K (X a topological space) that is proper and bounded below attains its greatest lower 
bound, which is therefore a minimum and hence a critical point. To see this, let ri be 
the greatest lower bound of f{X), and let r 2 > ri. Then K := f~^{[ri,r 2 ]) is non-empty 
and compact (as / is proper). Therefore f{K) is a compact subset of R so must be equal 
to [ri,r 2 ], and hence the greatest lower bound ri is attained by /, as claimed. 

Applying this to the restriction of fo to (/'o^(O) C Ox A^, for a = p+e sufficiently small 
that O CU, shows that indeed this restriction of fo has a critical point as required. □ 

Remark 4.7. There are two differences with the persistence theorem in [13]. The first 
is that in [13] it is only assumed that the relative equilibrium is ‘extremal’, rather than 
formally stable. In practice, formal stability is the most straightforward way of testing 
for extremality, although it is of course a stronger condition. An example of an extremal 
relative equilibrium which is not formally stable is the famous Thompson heptagon in 
the planar point vortex problem, see the proof in [14] of its stability (Corollary 14.7 and 
the calculation in Sec. 14.3.5.2). 

The second difference is more subtle. The conclusion in [13] is that relative equilibria 
exist for all p' in a neighbourhood of 3{z) in 3{Uq) (where, as in the proof above, Uq 
is the G-invariant neighbourhood on which the MGS normal form is valid). Here on 
the other hand, existence is only guaranteed for p' G 3(U), the image under J in 3(Uo) 
of some (small) G-invariant open set U. This image 3{U) may not in general be a full 
neighbourhood of p, = J{z) in J(17o)- On the other hand, if G^ is compact, then in [13] 
it is proved that J is G-open, meaning that the image of U will be open in J(17o) and 
the results will coincide (in this respect). 

5. Bifurcations 

We now consider the problem of bifurcations of branches of relative equilibria. Starting 
with a parametrized branch of relative equilibria of the same symplectic type, we will give 
sufficient conditions for the existence of bifurcating branches. At the end of this section we 
study the possibility of bifurcations from persisting branches consisting of formally stable 
relative equilibria. Recall that we assume G acts properly on the symplectic manifold 
V (the phase space) with momentum map J : 7^ —J g*, and this is equivariant with 
respect to a (possibly modified) coadjoint action. Moreover, h : V ^ M. is a G-invariant 
Hamiltonian function. 

5.1. Bifurcations from branches of constant symplectic type. 

The following result gives sufficient conditions for the occurrence of bifurcations from a 
parametrized branch of relative equilibria of the same symplectic type. 

Theorem 5.1. Let z €V be a relative equilibrium with momentum 3{z) = p and velocity 
f € Qf_i written ^ + V G ni© 0 z according to (2.1). For a subgroup K < (Gz)^ 

and a vector subspace W C m* x let z : W -P- V, be a parametrized branch of 
relative equilibria of the same symplectic type, satisfying ^(0,0) = z, 3{z{p,r]')) = p + p, 
Gj(z{p,rj')) =Gp, and Gz{p^r)') =Gz- Let £_{p,ri') be a family of velocities for points of this 
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branch chosen such that rj') = r] + rj' (this is always possible). Now suppose that 

there exists a subgroup L < K satisfying 

(i) 02 is co-central in 0 ^ (see Lemma 4-2), 

(ii) Nq^{L) acts on TV := with cohomogeneity one, and 

(iii) the eigenvalue a{p,r]'), with {p,ri') € W, o/ crosses 0 at 0 € W. 
Then, for every v € Af close enough to the origin, there is a relative equilibrium Ty near 
z with velocity G 0 p close to ^ and not in the G-orbit of any point of the original 
branch z{p,r]'). The stabilizer ofzy is 

Gz„ = (Gz)y > L. 


Proof. In the local model, the given branch z{p, rj') can be written as 

^9,1) = [9{p,'n'),P{p,'n'),v{p,'rj')], 

with [(/(0,0), ;o(0,0), o(0, 0)] = [e,0, 0], for (p, r?') G m* x 02 ^. Since relative equilibria 
come in group orbits, we can choose in our analysis g{p, rj') = e without loss of generality. 

Since points in the branch satisfy 3{z{p, rj')) = /i + p{p, rj') + 3ff(v{p, rj')) = p-\- p ii 
follows that 

P(P: P')= P and Gp+p = Gp 

Since z{p,r]') and z are of the same symplectic type, we can assume that Ggjp .,^ 1 '^ = Gz 
and in particular 

iP,vip,rj')) e X 

So actually W C X 0 ^^ C X 0 j,^. Let (pp^rj' : iV —t fV be the IL-parametrized 

family of diffeomorphisms of N given by 

Ppw'i'^) =V-\-v{p,Tj'). 

These maps are obviously iG-equivariant. Notice also that 3N{Pp,rj'{v)) = 3 n{v) since 
v{p, rj') G 

Notice that dhj(pp,/)(z(p, 77')) = 0 if and only if d/i^+p/(p, u(p, 77')) = 0 which, putting 
h'.pj^p,{p,v) = h{p,pp^r,'{v)) - {(fip^n'iv)), is in turn equivalent to 

(p, 0) = 0. 

Let us call f{p',p,v) = iip_|_p/(p, u). The function / G G°°(IL x N) is by construction 
(G 2 )p-invariant. We have according to Lemma 3.5 (using h' as the function h in the 
statement of the lemma) that, for every (p, 77') G W, 

Djvf(rj',p,0) = 0, 


and that for any L < K, the eigenvalues (and their multiplicities) of 

D%f{v',p,0)\^, 

and 


2 L I 

z{p,v’)"'iip,v') lAr-t- 


coincide. In particular, we also have 

keTD%f{0,0,0)\j^L = ker D%hp{0,0)\^L = kerd^/i^l^i. 
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We have that / satisfies the hypotheses of Lemma 2.9 with K = {Gz)jj. Therefore we 
can conclude that for each v G ker 0) |jYt\{0} there is a pair {p{v), rj'{v)) G W 

satisfying 

DNhr,+r,'(v){p{v),v) = 0 . 

We have shown that the point [e,p(f),u] satisfies (2.17). We now show that it also 
satisfies (2.16). The argument is again very similar to the one used in the proof of 
Theorem 4.3. Using the splitting 0 ^ = m* © 0 ^ we have that p{v) + JAr(j') G 0 *. In 
fact, since {p{v),v) G x N^, we can say that p{v) + Jiv(i') G 0 *'^. Calling 7 = 

Djn*h{p{v), v)+p + p', if we show that 7 G 0 ^ then (2.17) will follow from the hypothesis 
(i). It is clear that Dm*h{p{v),v) + ?7 + ? 7 'Gm© 02 = 0 ^, so it only remains to show 
that this element is invariant by L. This is true since, on the one hand p + p' € 2z 
on the other hand, using the Principle of Symmetric Criticality, an argument similar to 
the one used in Theorem 4.3 shows that Dm*h[p{v),v) G since p{v) G C m*^. 

Therefore for each v G kerj\/’\{0} close enough to the origin we have that z(v) = 
(/?([e, p(?;), ?;]) is a relative equilibrium with velocity Dm»h(p{v), v) + p + p' € Qfj,. In addi¬ 
tion, since = (Gj;)p(„)n(Gz)„ and we found that W G X 0 ^ then (Gj;)p(„) = Gz 
and therefore G^(„) = {Gz)v Since v G we also have L < G^yy □ 

5.2. Bifurcations from a branch of formally stable points of constant symplec- 
tic type. 

As a consequence of Theorem 5.1 we now show that branches of formally stable rela¬ 
tive equilibria of the same symplectic type with 2 -dimensional fixed point subspaces and 
satisfying the usual group-theoretic conditions used in this article, typically have bifur¬ 
cating solutions at each point of the branch. This is a phenomenon due exclusively to 
the existence of continuous isotropy and cannot happen for branches of formally stable 
relative equilibria with discrete stabilizers, for which the existing persisting branch of 
relative equilibria of the same symplectic type must persist without bifurcation. 

Theorem 5.2. Assume the hypotheses of Theorem 5.1 for a branch z{p,p') consisting 
of formally stable points, where conditions (ii), (in) are replaced with 
(ii) dimlV^ = 2 , and 
(hi is definite. 

Then, generically, for every v G kerd^/ij|^i\{0} close enough to the origin, there is a 
relative equilibrium 

Zy 

near z with velocity of the form G 0 p close to ^ and not in the G-orbit of any point of 
the original branch z{p,p'). 

Proof. Using the function / G G°^{W x N) introduced in the proof of Theorem 5.1 we 
have that in the local model we can assume that z{p, p') = [e, p, 0 ] and that each of these 
points is a relative equilibrium with velocity +77 + 77 '. Furthermore, the signatures 
of U>^/( 77 ',p, 0 )|j^i and ^^{pw')\n^ are the same. Since is definite then 

so is Dljf{r]',p,fi)\j^L in a neighbourhood of (0, 0) G W . We therefore have a map 
■0 : lU —J Sym(2) from a neighbourhood of (0,0) in W to the space of symmetric 2x2 
matrices and its image lies in the set of definite matrices. 
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If we give coordinates {x, y, z} on Sym(2), where x, y are the diagonal elements and z is 
off-diagonal, the set of definite matrices is given by the two connected components of the 
set Sym(2, -I-) defined by xy — z^ > 0. These connected components are the set of positive 
definite and negative definite matrices respectively. The set Sym(2, —) of symmetric 2x2 
matrices with index 1 is defined hy xy—z"^ < 0. These two regions are separated by the set 
Sym(2, 0) of matrices with at least one zero eigenvalue, given by xy — z^ = 0. The origin 
in represents the zero matrix. Therefore any point in Sym(2,0) other than the origin 
represents a matrix with only one zero eigenvalue. Suppose that there is an element fj £ K 
such that I?^J^(0) I is not proportional to Dj^f{r]', p, 0) | This implies that we can 
find some t G M satisfying £)^/(0,0)|^i G Sym(2,0)\{(0, 0,0)}. Notice 
that this last matrix is equal to D^/(0,0)|^i where f{ri',p,v) = hfj+tfj+rj'{p,v). Notice 
also that we can rewrite our starting branch of relative equilibria as z(p, p') = z{p, tfj+p') 
but now this produces the map 'ijj : W ^ Sym(2) given by '0(p,r?') = Dj^f{p',p,0)\j^L 
which satisfies i/’(0, 0) G Sym(2,0)\{(0,0,0)} and generically will cross this surface. In 
other words, substituting the initial velocity of z, the element by the also admissible 
velocity + p + tp, we can force all the conditions of Theorem 5.1 to be satisfied 
generically, and the existence of the stated bifurcating branch follows. □ 

6 . Example. The sleeping Lagrange top 

This classic example will allow us to illustrate our results on the persistence, stability 
and bifurcations of relative equilibria. We will briefly summarize the construction of this 
system, referring to [15] for more details, since we will be using the same notations and 
conventions. 

The Lagrange top is a mechanical system consisting in a rigid body with an axis of 
symmetry and a fixed point in the presence of an homogeneous vertical gravitational 
field. Relative equilibria are steady rotations around the vertical axis and the axis of 
symmetry. Mathematically, this is the symmetric Hamiltonian system 

(T*SO(3),a;„T2,J,h), 

where ujc is the canonical symplectic form on the cotangent bundle T*SO(3) and J : 
T*SO(3) —> its canonical momentum map. We will use the right trivialization of 
T*SO(3), and will identify so(3) with and with R^ in the usual ways. The Euclidean 
structures on these two vector spaces will allow us to identify them also with the duals 
of their Lie algebras. In this setup, a point in T*SO(3) is represented by a pair (A, tt), 
where A G SO(3), viewed as a matrix group, and tt G R^. The several relevant elements 
for the analysis are: 

• The action x T*SO(3) —>■ T*SO(3), given by 

((01,02), (A, t)) ^ (6.1) 

• The Hamiltonian 

h(A, tt) = mgles ■ Aea -|- — tt ■ 

where I = diag(Ji, Ii, Is) and Ia = AIA“^. Here Ji, Is are the principal moments 
of inertia of the body, I 3 corresponding to the axis of symmetry. 
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• The momentum map 

J(A, tt) = (tt • 03 , —TT • Aea) £ ~ 

The sleeping Lagrange top is a family of relative equilibria of this Hamiltonian system, 
consisting in the top rotating around the vertical position, with the gravity and symmetry 
axes aligned. This corresponds to the points in phase space of the form • z, where 
z = (A, tt) = (/, A/aea), I is the identity 3x3 matrix and A is any non-zero real number, 
corresponding to the angular velocity of the rotation. It is easy to see from (6.1) that 
the stabilizer of any sleeping Lagrange top is 

Gz = {(01,02) : 9i = 02}c^S\ 

The Lie algebra Qz can be identified with R as 

0z = {(»7,?7) £ : ry € R}. 

Therefore this relative equilibrium has continuous isotropy and its velocity is not uniquely 
determined. It is easy to see that a possible orthogonal velocity is given by = 
5 (A, —A) £ R^, corresponding to the choice 

in = {( 7 , — 7 ) £ R^ : 7 £ R}. ( 6 . 2 ) 

It follows that the set of all possible velocities for this relative equilibrium is (77 -I- 
A/2, r] — A/2), for any real number 77. Its momentum value is 

77 = (A/ 3 ,-A/ 3 ). (6.3) 

Since the symmetry group is Abelian, G^ = G = T^, and dimm* = 1. 

It follows that with respect to the basis of N given by 

■Cl = (ei, 0 ), V 2 = (e 2 , 0 ), 773 = ( 0 , ei), V 4 = ( 6 , 02 ) 

the symplectic normal space N a.t z = (/, A/ 303 ) can be identified with R^, and 

d20(r,+A/2.r)-A/2) = 

where 

A 2 J 

A = -mgl - ( 2/3 - /i) - 77 A/ 3 , 

All 

Stability. We will fix A 7 ^ 0 and study the nonlinear stability of z = (/, A/ 3 ). In order 
to apply Theorem 3.6 we need to find an admissible velocity ^ for z such that is 

definite. This velocity will be of the form ^ = (77 -I- A/2 ,77 — A/2). 

The eigenvalues of dlh(^^+x/ 2 ,r,-\/ 2 ) |jv are given by 

^± = 1 ((Al + G) ± y/(A + Cy-4(AC-B^)^ , 


(AO OB 
0 A -B 0 

0 -B C 0 

\ B 0 0 C ) 


(6.4) 


B = A(273-7.)-,, C = j-. 


(6.5) 
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each with multiplicity 2. It follows that this matrix is definite provided AC — > 0. 

Since ij is arbitrary and A 7 ^ 0, we will make rj = kX for arbitrary fc G R. Then we have 

AC - = ^ (A^ {h{Ak + 2) - /i(4fc2 + 4fc + 1)) - Aglm) . 

47 1 

Then, the relative equilibrium will be stable if we can find some fc S M for which A^ > 
/ 3 ( 4 fc+ 2 )-Zi( 4 fc^+ 4 fc+i) • clear that this is always possible if A^ is sufficiently large, so 
we are interested in the element of the family of sleeping Lagrange tops with minimum 
|A| yet nonlinearly stable. For that, we find a minimum of the above function of k, which 
happens when k = and then, the relative equilibrium z = (/, A /3 63 ) is nonlinear ly 

stable provided 

2 ^ Aglmh 

^ > —72— 

which is the well known classical stability condition (fast-top condition). 


Persistence. It is straightforward to show that the linear action of Gz ^ on N is 
given by 


/ COS 0 —sin 6 * 0 0 \ 

sin 6 cos 0 0 0 

0 0 cos 9 — sin 9 

\ 0 0 sin 0 cos 0 / 


( 6 . 6 ) 


It follows that = {0}, whence is trivially non-degenerate. Note also 

that is a torus. Recall that (T*SO( 3 ))( 5 i) is the set of points of the form 

(/, XI^Bs) with A 7 / 0, that is, the set of all Lagrange tops which are not at rest. Taking 
into account expression (6.3), it follows from Theorem 4.4 that near a sleeping Lagrange 
top with phase space point (/, Ao/sCs) there is e > 0 such that for every r £ (—e, e) 
there is a T^-orbit of relative equilibria of orbit type (S*^) that forms a persisting branch 
of nonlinearly stable relative equilibria of the form (/, (Aq -f r)/ 3 e 3 ). Furthermore this 
branch is a local symplectic submanifold of T*SO(3). 

Note that in this particular case, the persistence result already follows from Remark 
4.3 in [5], although in that reference the conditions are stronger in general, requiring the 
relative equilibrium to be non-degenerate. By the uniqueness argument of Theorem 4.4, 
it follows that this persisting branch corresponds exactly to a neighbourhood of A in the 
family of sleeping Lagrange tops, which exists for any non-zero A, as can be shown by a 
direct computation. Note that the case A = 0 is not contained in the persisting branch 
of sleeping Lagrange tops since in this case = T^, and then the symplectic type of 
this point is different from the elements of the branch. The case A = 0 corresponds to 
the unstable equilibrium of the top in the upright position. 


Bifurcations. We now apply Theorem 5.1 to study the possible bifurcations of rela¬ 
tive equilibria from the family of sleeping Lagrange tops. We will start by a relative 
equilibrium in the family of sleeping Lagrange tops of the form z = (J, Xohea) with mo¬ 
mentum g = (A0/3,—A0/3) and orthogonal velocity = (Ao/2 ,—Ao/2 ) £ m according 
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to the choice (6.2) . Notice that since Gz is Abelian, (Gz)t, = Gz — . Let us choose 

L = K = {e} and W = m* x {0} Cm* x We can identify 

^ = {(p,-p) : P e K} 

and define a local branch of relative equilibria z(p) = (I, (Aq + p/l 3 )l 3 es). Notice that 
z(0) = z and 3{z{p)) = p,+ {p, —p). It is also clear that Gj(j(p)) = p+{p, —p) = Gfj, = G 
since G is Abelian, and also that Gz(p) = Gz — S^. We choose the family of velocities 
C(p) = ((^0 + ^)/2, —(Ao + r)/2) + ( 77 , 77 ) for a fixed, although arbitrary, element 77 G K. 

It is clear, again from the Abelian character of G, that condition (i) in Theorem 5.1 
is trivially satisfied. In order to check condition {Hi), notice that = N. According 
to (6.5) we see that cr+{p) > 0 for every value of p but 

sign(tT_(p)) = sign(AG — B^). 


Computing the value of this last expression along the branch z{p) we get 
AhlU^C-B^) = p'^{2l3-h)+p{{2r] + Xo)l3{Ii-l3)-2I^Xo) 

- ((277 + Ao)^/i /3 - 2(277 + Ao)/|Ao + ‘imgll^) . 


It is straightforward to check that 77 can be chosen such that this expression vanishes 
at p = 0 if one can solve the independent term for 77 . This is possible precisely when 

Amglli 


K> 


^3 


that is, if the local branch z{p) is centred at a point in the formally stable range. In order 
to check that the eigenvalue ct_ (p) actually changes sign at p = 0 we have to guarantee 
that the coefficient of the linear term is different from zero. But this is true since if both 
the coefficients of the linear and independent term vanish simultaneously we would have 


2 I 3 X, 


3^0 


h — h 


hll - 2 


2 fyAo 

h —13 


/|Ao + Arngll^ = 0 


and the above expression is always positive so cr{p) changes sign at p = 0 and condition 
(Hi) in Theorem 5.1 is satisfied. 

Thus it remains to study condition {ii). For that, notice that since L = {e} and 
Nsi{e) ~ SO(2) we have to check if 




with 77 chosen in a way that the crossing condition (Hi) is satisfied, is two-dimensional 
and invariant under the action of Gz — SO(2) given by ( 6 . 6 ). It is easy to see from (6.4) 
that for A = Aq, 


kerdj/i(^+A/ 2 .r,-A/ 2 )|^ = span( ( a, 0 , 0 ,--^a 


0 , 6 , ^ 6,0 


with respect to the basis {ui, 712 , 773 , 774 } of N. It is clear from ( 6 . 6 ) that this space 
is a S'^-module equivariantly isomorphic to equipped with the standard action of 

S^. Therefore (ii) is also satisfied. It follows from Theorem 5.1 that for every v G 
kerd^/ 7 (^_i_A/ 2 ,r)-A/ 2 ) Ijv near the origin, there is a relative equilibrium for the Lagrange 
top system not contained in the branch of sleeping tops. Since the action of on 





HAMILTONIAN RELATIVE EQUILIBRIA 


33 


this space is free outside the origin, these relative equilibria have trivial isotropy and 
correspond to processing tops, where the symmetry and gravity axes are not aligned. 

7. Example. 2 point vortices on the sphere 

Consider the system of 2 point vortices on the sphere [4, 7], which is a simple system 
where every motion is a relative equilibrium. The phase space is P = x S^\ A, where 
A is the diagonal, with symplectic form a; = Tiwi +T 2 UJ 2 , where uij is the natural SO(3)- 
invariant symplectic form on the copy of the sphere and Tj G R are the vorticities of 
the two points (taken to be non-zero). The equation of motion is 

„ X2 X Xi . Xi X X2 

xi=r2--, X2=ri- -. 

1 — Xi • X2 1 — Xi ■ X2 

This is a Hamiltonian system, with Hamiltonian 

/i(xi,x2) = -rir2iog||xi -X 2 IP 

where ||xi — X 2 II is the Euclidean distance between the points (different versions of this 
differ by factors of 2?!, but this can always be compensated by rescaling time). There 
are two possible G-Hamiltonian systems depending on whether Ti or r 2 are different or 
equal. If Ti ^ r 2 then the group SO(3) acts on V by the diagonal action. We think of 
each copy of as embedded as the unit sphere in This action is Hamiltonian with 
momentum map given by, after identifying so(3)* with R^ as usual, 

J(X1,X2) = TiXi -b r2X2. (7.1) 

IfTi = r 2 then the system supports a Hamiltonian action of the direct product SO(3)xZ2 
where SO(3) acts as in the previous case and ZJ is the reflection group generated by 
t(xi,X 2 ) = (x 2 ,xi). The momentum map for this action is also given by the expression 
(7.1). 

The action of SO(3) (respectively SO(3) x ZJ) on phase space is locally free except 
where xi = —X 2 (antipodal points), in which case the isotropy is SO(2) (respectively 
0(2), which is generated by the group S0(2) of rotations around Xi and an element 
{Pn,T) £ S0(3) X Z 2 . Here Rn denotes a rotation of angle tt around a vector n per¬ 
pendicular to xi. Since, for a given S0(2), the subspace consists of just 

two points, they are necessarily equilibria for any invariant Hamiltonian. Non-antipodal 
pairs have stabilizers {e} or Z 2 (n) respectively, where Z 2 (n) is the subgroup generated 
by the element (i?n, t) G S0(3) x Z^. 

With 2 point vortices, as already mentioned, every solution is a relative equilibrium, 
and indeed given any xi,X 2 which are not antipodal, the angular velocity is 

^ “ 2sm^{e/2)^ 

where = J(xi,X 2 ) and 9 is the angle subtended by the two points. With respect 
to antipodal points, the situation changes depending on the two possibilities for the 
vorticities. 

1. Case Ti ^ r 2 . In the limit of an antipodal pair, so as X 2 ——xi, one has ^ t = 
i(ri — r 2 )xi ^ 0. This limiting ^ is of course an element of Qz for z = (xi,—xi), in 
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order that z be an equilibrium point. In this case N is 2-dimensional. Specifically, take 
z = (es,—es) S S'^ X and using coordinates {xi, X2,yi,y2) the tangent space to the 
group orbit is spanned by (1, —1,0, 0) and (0, 0,1, —1). Since dimg^ = 1 the rank of T^J 
is also 2 , and its kernel is spanned by the tangent vectors 

Vi = (r2, -bi,0,0), V2 = (0,0,r2, -bi). (7.2) 

This is therefore also the symplectic normal space N when bi ^ b 2 . On this space, and 
with this basis, 

dz/lU = ^rib 2 (bi - b 2 ) 2 j^, d 2 j|^ = bib 2 (bi - b 2 )/N, 

where In is the identity matrix on N. Thus, for ry S ~ M, 

dlh^lN = ^rib 2 (bi - b 2 )(bi - b 2 - 2r,)lN. 

Thus, assuming bi ^ b 2 , this is degenerate precisely at 77 = ^(bi — b 2 ); note that 
this is precisely the limiting velocity mentioned above. In particular, note that 
yi = (bi — b 2 )e 3 ^ 0 and then = Gz = SO(2) and m* = 0. 

With reference to Theorem 4.5, the relative equilibrium z is an equilibrium, so = 0, 
and let 77 = ^(bi — b 2 ). We take L = {e} and K = Gz = SO(2). Therefore = N and 
= 0 since G^ is Abelian. We have 

dzS-l-V Ijv = rib2(b2 - bi)77'/Ar, 

which has a unique eigenvalue = bib 2 (b 2 — bi) 77 ' that changes sign at 77 ' = 0 , 

satisfying hypothesis (*). At 77 ' = 0, kerd^/i^j^ = N which is 2-dimensional, and from 
(7.2) is a SO(2)-module isomorphic to equipped with the standard circle action, 
satisfying hypothesis {ii). Since the action of Gz is free outside the origin, it follows from 
Theorem 4.5 that for every element of A^\{0} there is a relative equilibrium near z with 
trivial isotropy. That is, a rotating configuration of non-antipodal pairs. 

On the other hand the content of Theorem 4.3 is empty in this case, since m* = 0 
and therefore we only get the same equilibrium point, with different velocities t]' € Qz- 
On the third hand. Theorem 4.4 is not applicable here since P(so( 2 )) consists of a single 
group orbit. 

2. Case bi = b 2 . In this case, noting that for z = (xi, —Xi) we have y, = J(z) = 0 and 
we find that G^ = SO(3) x ZJ and Gz = 0(2). Therefore the kernel of T^J coincides 
with the tangent space to the group orbit at z, which implies that N is trivial. However 
m* = xy*- and is therefore 2-dimensional. Since is trivially non-degenerate for 

any 77 ', Theorem 4.5 is not applicable. Also, since in this case Gz = 0(2), then 
consists of a single group orbit and Theorem 4.4 is not applicable either. 

With respect to Theorem 4.3, let us choose K = Z 2 (n). Then we have m*^ = Mn 
which is Abelian, fulfilling condition ( 7 ). We also have = 0. Since is trivial, 
condition (ii) is automatically satisfied too. It follows from Theorem 4.3 that there is 
a smooth branch of (G-orbits of) relative equilibria near z with stabilizers containing 
Z 2 (n). Noticing that for bi = b 2 the stabilizer of non-antipodal pairs is precisely Z 2 (n), 
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we find that the application of the theorem produces the branch of rotating non-antipodal 
relative equilibria near antipodal equilibria. 

Since every motion is a relative equilibrium, the reduced dynamics for this system is 
trivial so stability holds trivially. Moreover, the theorems of Sections do not apply as 
there are no non-trivial branches of relative equilibria from which to bifurcate. 
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